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In this paper, the principles of the general relativity are used to formulate quantum wave equations 
for spin-0 and spin- 1/2 particles. More specifically, the equations are worked in a Schwarzschild- 
like metric. As a test, the hydrogen atom spectrum is calculated. A comparison of the calculated 
spectrum with the numerical data of the deuterium energy levels shows a significant improvement of 
the accord, and the deviations are almost five times smaller then the ones obtained with the Dirac 
theory. The implications of the theory considering the strong interactions are also discussed. 



I. INTRODUCTION 



• The general theory of the relativity, proposed in 1916 by Einstein, was one of the major scientific discoveries of last 
century. Besides providing very accurate theoretical results, it was a great advance in the understanding the Nature, 
£^ ■ dealing with the structure of the space-time. 

1} ' A question of interest is how the quantum theory can be affected by the space-time. Dirac, formulated his theory 
, based in the flat space-time of the special theory of the relativity, and with this formulation, the spin and the 
antiparticles appeared naturally into the theory. 

About curved spaces, many authors, as for example 0, 0, Q proposed methods to quantize the gravity, where 



& . still there are many difficulties and opened questions to be understood. 



In this work, a different point of view is proposed. Here, instead of trying to quantize the gravity, the effects of the 
metric in the subatomic world will be studied. For this purpose, the basic idea is to describe a particle in a region with 
a potential that affects the metric of the space-time. We are not interested in gravitational effects, as in 

0- 

where the 

effect of gravitational forces in the hydrogen atom spectrum has been included. So, the gravitational potential will 
be turned off and only the other interactions (strong, electromagnetic) will be considered. Observing that the masses 
of the particles are very small, and the small value of the gravitational coupling, when compared with the electric 
or strong ones, one can say that it is an excellent approximation. Inside this space-time, curved by the interaction, 
according to the general covariance, quantum wave equations will be proposed. Then some simple applications will 
Q\ , be made, in order to verify the predictions of the theory. 

This paper will show the following contents: In Sect. II the operators in the Schwarzschild metric will be calculated, 
in Sect. Ill, a brief review of the dynamics will be made, in Sects. IV and V the quantum wave equations will be 
proposed. In Sect. VI we will apply the theory to the hydrogen atom, calculating its energy spectrum and in Sect. 
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C/3 | VII, the strong interactions will be considered. In Sect. VIII, the conclusions will be presented. 

>V II. THE METRIC 

Oh 



In this section we will calculate the operators (E, p, p 2 ) needed in order to write the wave equations, using the 
general relativity principles. As a first step, a system with spherical symmetry will be considered, but the basic ideas 
can be generalized to systems with arbitrary metrics. 

We will consider a particle inside a field, that may be described by a potential function V. The source of the field 
(a mass for a gravitational field or a charge for an electromagnetic field) will have some distribution, described by a 
tensor T^ v ^ 0, in a certain space region. 

Outside of the source distribution, on the empty space (where T^ v ^ 0) if we consider a system that presents 
spherical symmetry, with a central potential V(r), the space-time may be described by a Schwarzschild- like metric 
0,0, 

ds 2 = £ dr 2 -r 2 (d0 2 + sin 2 6» d<p 2 ) - ^dr 2 , (1) 



where £(r) is determined by the interaction potential V(r), and is a function only of r, for a time independent 



2 



interaction. £(r) will be studied in detail in Sec. III. As we can see in QJ, the metric tensor is diagonal 
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and can be written in the form 



where 
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Using these definitions, we can calculate the operators 

p = —ih\7 . 

According to the above expressions, the momentum operator may be defined as 

6 8 



p = —ih 



. /-a 9 8 



' dr r 89 r sin 9 

that in a region with V = (£ = 1) is the usual momentum operator in spherical coordinates 
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p = —ih 



„ a e d 

dr r 89 ' r sin 9 dip 



The energy operator is defined as 
that for V = 0, 

The Laplacian is calculated using 

V 2 = (h^hs) 
where hi are given in J5J), so 

\pf = -h 2 



E = ihW 



ih d ih d 
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(2) 



(3) 



(4) 
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If one defines the momentum components as pi — —ihVi, one can observe that they are not good operators, these 
operators does not commute and are not even Hermitians, so the definition [T(i| 



(13) 
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where D 



-g, with g = det(g fll ,), will be used. With this definition, 



and the commutation rules are 



The Laplacian can be expressed as 



V 2 = 



that for weak potentials is just 
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(14) 
(15) 
(16) 
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With the operators calculated in this section, one can obtain the relativistic quantum wave equations. If another 
symmetry is important (as axial symmetry, for example ) , the operators can be obtained in a similar way in the given 
metric. 



III. SCHWARZSCHILD DYNAMICS 



In order to obtain the wave equations, two expressions are needed: the energy and £(r). Then, also with the 
function of setting the notation, it is useful to make a brief review of the dynamics, in the Schwarzschild metric, and 
to show how the quantities of interest can be expressed. 

From the proper time is 



dr 



<ds 2 



ft 



drJZ-^ + r 2 (3 2 



with 



7s = 



£-f + r 2 2 



where (3 r and fit are the the radial and transverse parts of f3, respectively. They are defined as 

B = - 

1 dr ' 

(3r = 



dr 

dr. 



The principle of least action states that 

S = [ Ldt 



-ITIqC 



ds 
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(21) 
(22) 

(23) 
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where mo is the rest mass of the particle. Then, the Lagrangian can be expressed as 



L = -m c 2 ] Jt;- y+r 2 /? 2 =-m Q c 2 / ls . (27) 

The momentum four-vector is defined as 

p» = Bop = E Qls {lJ) = (p°,p) (28) 
^ = (£/,-rV,-ry,-r 2 sin 2 0p 3 ) , (29) 

where 

E = m c 2 . (30) 

The equivalence principle provides the relation 

dr 

that gives 



^- = -T%^^ , (31) 



ls d ^ r = E TIM V = E [T^p + = . (32) 



Po = E = ^ S E = = , (33) 



So, the energy defined as 



t 

is a constant of motion. The other constant is L z = ps/c. 
In the rest frame of the particle 

p 0f j% = -^o = -m 2 c 4 , (34) 

that is a Lorentz invariant, so 

p^pV-y = -mgc 4 (35) 

and then, the energy relation is 

£ 2 



p 2 c 2 + m 2 c 4 , (36) 



or 



— = ^p2 C 2 + m 2 c± . (37) 

The expressions (|36l) and (j37(l will be used to construct the Hamiltonian operators. 
The term £ is a function of r, and can be determined if we observe (|37|l 

E(P = 0) = £ £ 1/2 = + V , (38) 

that means that in the rest frame of the particle, the energy is due to the sum of its rest mass with the potential. 
Then 

e /2 = 1 + ^ = 1 + ^ ■ (39) 
Eq m c^ 

Comparing with the standard definition of the Schwarzschild mass 

2 m s 2V V 2 

e = i - = i + — z + — , (40) 
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it is possible to make the identification 



r ( 2V V< > 
ms = -2{^ + ^) ' (41) 



that for general potentials may be a function of r. 
For weak potentials, 



2V 

£-1 + - . (42) 



moc 2 



and 



Vr 

m s ? (43) 



TOqC 2 



that are the usual expressions of general relativity, 



/ GM\ 2 2GM , 

= i a ~ i r • 44 

V r cr ) r cr 

IV. SPIN-0 PARTICLES WAVE EQUATION 

With the knowledge of the energy 136fl and £(r) (|39|l . it is possible to formulate wave equations in the given metric. 
The simplest case is to obtain the equation for spin-0 particles. For this purpose, the procedure to be followed is the 
same one that is used to determine the Klein-Gordon equation, that is based in an operator for E 2 . Using the relation 



i? 2 _ 

and H36f) the quantum wave equation, based on general relativity, for spin-0 particles is 

h 2 d 2 ^ 



(45) 



e dt 2 



= -h 2 c 2 ^ 2 ^ + mf,c 4 ^ , (46) 



with V 2 defined in (JTTJ. 

This equation can be separated in the standard way, yielding 

y(r,6,^t)=u(r)Yr(8,<l))e- lEt / h , (47) 

where Y l m (9, <p) are the spherical harmonics. The radial equation is then 

' -*yflsz ) u +[ tttto - -4 —r 2 - ) u = . ( 48 ) 



r 2 dr y T v<, drJ " ' \h 2 c 2 £ 2 fi 2 
and can be solved for a given interaction potential V(r), that determines £(r). 

V. SPIN-1/2 PARTICLES WAVE EQUATIONS 

The next step, is to obtain the analog of the Dirac equation, for spin-1/2 particles. It can be made in the same way 
that it was made by Dirac using an Hamiltonian with the a and /3 Dirac matrices, instead of an energy operator 
with an square root (|37|l . Then we have 

= i-ihc a.V + (3m c 2 ) * , (49) 
5 ot V / 

and if we square the operators in both sides of the equation, we obtain the wave equation l|46|). that proofs that the 
procedure used by Dirac, is also valid in this case. 
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Separating the time dependent part 

T(t) = A e - lEt/h , (50) 



we will have the spatial equation, 



E 

-ihc a.V + /3m c 2 -= ) t/j(r) = . (51) 



Observing the relation 



a.V = —ol.t 



r 



-1/2. 9 , „r d t a* d 



dr \ r 86 r sin 



(52) 



and the angular part of the operator of (|52|l . one concludes that the angular part of V'(^) can be described in terms 
of the two component spinors xt' 13 > 

Xl= ]T C(l,l/2,j^-m,m)Yr m (0,<t>)x m , (53) 

m=±l/2 

where C(/, 1/2, j; /i — m, m) is a Clebsh-Gordan coefficient, x m , a Pauli spinor and 

k = l for j = 1- 1/2 , 

k = -l-l for j = / + 1/2 , (54) 

that gives 

fc = ±(j + l/2) . (55) 
Then, the wave function is expected to have the structure 

^=(1&I ) , (56) 



with the F and G functions obeying 



dF 




F 

hfe - = 
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dr 




dG 
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-(- 


dr 





mojF , (57) 

that are equations very similar to the ones obtained from the Dirac theory. In the following sections, some physical 
implications of this theory will be studied. 

VI. THE HYDROGEN ATOM 

In this section we will study the behavior of the theory, in a very well known system, the hydrogen atom. In the 
hydrogen atom, the electron is submitted to an electric central potential (obviously with Z = 1) 

. , aZ 

V(r) = , (58) 

r 

then the £ function becomes (where ep means electron-proton) 

£ ep = fl V . (59) 
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where m e is the electron mass. This function represents a space-time curved by the ep interaction, or how space-time 
is seen by the electron. At this point we can see an interesting feature of the theory: when general relativity is used 
to study the gravitation, the radius where the metric breaks (£ = 0) , that is the Schwarzchild radius, r», is always 
negligible, as for example r =2.95 Km for the sun. But in the hydrogen atom (with the parameters of [15J) 

= 2.818 fm , (60) 



m e c 2 



that is the classical radius of the electron, which is obviously not negligible. The estimated radius of the proton is 
about 0.9-1.0 fm, so, it is located inside the horizon of events. Then, the electric charge will be confined inside this 
region by a trapping surface, as defined in [l4|, and outside, only effects of the total charge can be probed by the 
electron, and no information about the inner structure can be obtained. 

Inserting £ ep from the expression (|59|l in the equations (|57JI we obtain the equations, but valid only for r > r s . 
Inside the horizon of events, the metric is not the same, the energy-momentum tensor determined by the charge 
and matter distributions must be considered. In this paper we shall study only the outside behavior. It is illustrative 
to study the approximation 

+ (l + k)— = — — + m e G 



dr v ' r V 1 + v / E o 

EV \ 
E - — +m e ) G , 
Eq J 

dG n ,n G ( E \ n 
h (1 - k)— = - ; m e ) F 

dr 1 'r \1 + V/Eq J 

- - ( E -^- m *) F ' ( 61 ) 

where we neglected the V/E and higher order terms. As we can see, in this theory, the Dirac theory is recovered 
only for E/Eq ~ 1 (lower moments), 

dF F 

— + (l + fc)- = {E-V + m e )G 
dr r 

dC G 

-T~ + (1 — k)— = — (E —V — m e ) F . (62) 
dr r 

We must remark that in the approximation of equation 1)61(1 , the metric divergence at r — r s (that is not a physical 
singularity, it appears from the choice of the coordinate system) is removed, then £ ~1 and r s does not exist. So, 
in this case, the wave functions must be controlled only at the physical singularity, at the origin. Considering the 
Frobenius method, the solutions are of the type 

N 



n=0 
N 

G = P S J2 h np n z~ P ■ (63) 



n=0 

where p = (3r, and after some manipulations one finds 



(3 = \/m 2 c 4 - E 2 (64) 

and 



-L- . (65) 
m l 

Observing the solutions (|63|l one can see that F(r ~ 0) and G(r ~ 0) are sums of terms of the type Ce~ p p s+m (r ~ 0), 
and then, ip{r — 0)=0. So, the effect of the approximation, is to remove the horizon of events, and to extend the 
solution to the region r < r s . Considering the exact solution of the equation, this behavior is valid only for r > r s , 
and an horizon of events exists at this surface with the proprieties described above. But one must remark that at 
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the atomic level, an approximation of the order of 2.8 fm is not so bad, in this region the wave functions are almost 
negligible and for practical purposes this approximation is reasonable. 

From these expressions, it is possible to calculate the energy spectrum. Using the standard methods ^jj, 
one obtains 



(m e c - E )a - — , 



(66) 



that has the solutions 



with 



and 



E 2 = m 2 c 4 



-1± v/l + 4 7 2 /a 2 



a 
?ic 



a = o(n) = n - (j + 1/2) + y/ (j + 1/2) 2 - 7 2 
Using only the positive root solutions of <|67|l the hydrogen atom spectrum is 



(67) 



(68) 



(69) 



E n =m e c 2 J . (70) 



Adopting the expansion (7 /a is small) 



^2 Jl ~4 6 8 10 

1 + 4' ^ 1 + 2-L - + - 1Q-L + 28^ + .... (71) 
a a a a° a a u 



the energy spectrum l|7U|) can be rewritten as 

E = meC2 (72) 

n 71 + 7 2 /a 2 - 7 4 /a 4 + 2 7 6 /a 6 - S-f/cfi + ... ' 

where we can find exp licitly the corrections of the energy levels, due to general relativistic effects, if compared with 
the standard |12| relativistic spectrum 

m e c 2 , , 

E n = , , (73) 

VI +7 2 /a 2 

that can be obtained from the Dirac equation or from the Sommerfcld model 

Considering now the spectrum obtained from the exact solution of 157fl . without the approximations made in (|61() 

^- mec2 V2-8^ + 4^Vi^ + 2 ' (?4) 

we will compare the theoretical results with the experimental data and with the ones obtained with the Dirac theory. 

Some numerical values are shown in Table I, the experimental results |l6[ for the differences between the energies 
E(n, l,j) and the ground-state energies Ex, for the hydrogen atom and deuterium, the corresponding values calculated 
with the Dirac theory l|73ll . and the results calculated in this work, with l|74l) . Observing the table, one can see that 
the accord of both theories with the hydrogen experimental data is very good, but the results from (|74|l are closer 
to the experimental data then the results from (|73|l . One must remark that spherical symmetry is not exact in the 
hydrogen atom as the proton mass is finite, but with a heavier nuclei, this symmetry is a better approximation, so it 
is interesting to observe the deuterium data. Comparing the results form the Dirac theory (|73|l , one can see a better 
accord, and the deviations from the data are of the order of 0.027%. Considering the spectrum l|74|) . the deviations 
are of the order of 0.005%, approximately five times smaller. 
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TABLE I: Experimental energy levels (eV) for the hydrogen atom, for the deuterium [l6|. and the theoretical ones, calculated 
with the Dirac theory l|73^ and with l|74[l . 





Hydrogen 


Deuterium 


Dirac 


Eq. lf74jl 




-13.59844 


-13.60214 


-13.60587 


-13.60298 


E(2, 0,1/2) -Ex 


10.19881 


10.20159 


10.20444 


10.20172 


E(3, 0,1/2) -Ei 


12.08750 


12.09079 


12.09413 


12.09127 


E(4, 0,1/2) -Ex 


12.74854 


12.75201 


12.75551 


12.75263 


E(5, 0,1/2) -Ei 


13.05450 


13.05806 


13.06164 


13.05875 



VII. STRONG INTERACTIONS 



In this section, the implications of the theory, when strong interactions are taken into account will be studied. The 
simplest system possible is the NN interaction. If, as a first approximation, only the long range part of the potential 
would be considered, it should be dominated by the one pion exchange contribution (Yukawa potential), 



2 e 



V(r)=g A (75) 

r 

where g 2 = 13.40 is the NN coupling constant and \i is the pion mass. As V(r) is a function only of r, if we locate 
one nucleon at the origin, we would have 



*N N 



i-g 2 ^-^) ■ (76) 



However, as it is well known, the NN potential is not central, and there are other contributions, such as the tensor 
part 0, 0, that arises from more complex processes (two pion exchange [2(| and others). In order to make 
some estimates, some symmetrical cases of the Reid ^3 potentials can be used, as they are phenomenological ones 
and can give a first idea of the Schwarzschild radius. The potentials are superpositions of Yukawa type terms, 

VQ-S) = -h (e- x + 39.633e- 3x ) /x , (77) 
VCD) = -h(e- x +4.939e- 2x + I54.7e- 6x ) /x , (78) 
VCS) S = - (he- x + 1650.6e- 4:E - 6484.2e~ 7x ) fx , 

(79) 

where x = 0.7 r fm~ , r is the relative radius and /i=10.463 MeV. Inserting these potentials in 139|) and equating it 
to 0, we will have 

r s CS) = 0.33 fm , (80) 
r s CD) = 0.44 fm , (81) 
r s CS) s - 0.33 fm . (82) 

Thus, the part of the source of the strong force that is inside the horizon of events will be submitted to the the 
trapping effect |T^|. that prevents the escape of any matter and radiation, what means radial collapse of the source 
of the strong forces. 

It must be noted that the potentials containing spin dependent terms and others were not considered, fact that 
would break the spherical symmetry. However, in a first approximation, these terms can be considered as corrections 
to the central potential. But even if we consider these terms, with another metric, as for example an axial symmetric 
metric, the trapping surface would still exists, confirming the present conclusions, and, only giving a more accurate 
estimate of the size of the confining region. Collapse is not an exclusive feature of spherical symmetric systems, as it 
was stated in |14). deviations from spherical symmetry cannot prevent space-time singularities from arising. 

As we can see, when strong interactions are considered, the horizon of events is located in a radius that is not 
negligible. The preceding example gives an estimate in the range of 0.3-0.5 fm. This fact suggests that the quark 
confinement may be understood from these results. To understand the mechanism, let us consider that the source 
of the strong force obeys some matter distribution. Each element of this matter distribution, (that may be a quark, 
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but in general, this assumption is not necessary) with mass too, suffers the action of an attractive strong force. Some 
models [2]]], |2^,|2^, HH consider central Coulombic potentials of the type (|58|l to describe the effective interaction 
to which the quarks are submitted inside an hadron. A good example is the Cornell model [2^], that with a linear 
plus Coulomb central potential 

V(r) = -- + ± , (83) 
r a z 

with the parameters a ~ 2.34 GeV -1 and a ~ 0.5, is able to describe the J/ijj and the T. 

Thinking in terms of constituent quarks it is possible to use the proposed theory to give a description of some 
hadrons. Table II shows an estimate of a (for a Coulomb potential) and of the constituent quark masses (to) in order 
to obtain the masses (M) of the nucleon and the J/xj) and T mesons. The experimental values of these masses are also 
shown. These constants define the value of r s , inside of which the quarks are expected to be confined. The values of 
r s = 0.83 fm for a nucleon and 0.05 — 0.11 fm for heavy mesons are very reasonable and show that for heavier quarks, 
the values of r s are smaller. Table II was constructed with the objective of giving an idea of the magnitude of the 
constants, but a detailed description of the observed hadrons must consider additional terms in the potential. 



TABLE II: Values of the masses M compared with the experimental ones for some systems. The calculations are made 
considering Coulomb potentials with the parameters a, m and r s . 





m(GeV) 


a 


r s (fm) 


M(GeV) 


M cxp (GeV) 


Nucleon(gqg) 


0.38 


1.60 


0.83 


0.938 


0.938 (proton) 


J/V(cc) 


1.79 


1.00 


0.11 


3.10 


3.10 


T{bb) 


5.50 


1.05 


0.05 


9.47 


9.46 



With these results, it is possible to calculate [2f|, |27 



where 



2± = 5 <az> =hl-2 5), (84) 
gv 3 3 



3J ' — 0.059 (85) 



/(|F(r)|2 + |G(r)P)dr 

where a nucleon composed of tree quarks with j z = 1/2 is considered. So, gyt/gy^l.47, what shows a 17% deviation 
from the experimental result that is 1.259. The magnetic moments of the proton and the neutron may also be 
calculated 

"> = il ^ L = 2 ' 82 

, n = -^^ = -1.88, (86) 

that are in good agreement with the experimental results /z p =2.79 and /i n =-1.91 [l5| . 

If the whole quark content of the hadron is located inside r s (now, spherical symmetry is a good choice), the 
classical description of such a system would predict the collapse of the whole matter in the singularity located at 
r = 0. However, as we are dealing with a quantum system, the uncertainty principle will prevent this collapse. 
Consequently, there are two opposite effects acting on the elements of matter, and the resulting effect will be radial 
oscillations. This effect, in a flat Minkowski space-time may be described by effective potentials of the form 

V e g ~ a + air + a 2 r 2 + ... . (87) 

There is no surprise why some authors pfij . p9|. explain the hadronic structure with models based on harmonic 
oscillator quark models or with linear potentials of the type Ar[3(J, (25|. Inside the horizon of events, V > Eq, so Eq. 
(|S"7Jl can be expanded as 
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where terms similar to the effective potential l|S7|) appears. So, considering a Coulombic potential (fact that is not 
strictly necessary, other kind of potentials may present similar proprieties) with the correct parameters, confinement 
effects may occur, fact that is widely used, with the addition of confining potentials in flat space-times. 

VIII. SUMMARY AND CONCLUSIONS 

In this paper, quantum wave equations, based on the general relativity, in a Schwarzschild-like metric, have been 
obtained. Investigating the hydrogen atom spectrum, the approximate expression resulting from the theory (|70J) is in 
accord with the experimental values, and shows a small improvement due to the general relativistic corrections, when 
compared with the standard relativistic spectrum l|73|l . Although, if the exact solution is considered, the corrections 
are not so small and (|74|l gives a significant improvement of the accord with the experimental data, specially with the 
deuteruim spectrum. 

An interesting feature of the theory, is that in the Schwarzschild metric, the horizon of events appears for r — r s , 
with a value that is not negligible, as it happens when the gravitational interaction is considered. When considering 
the strong interaction, r s shows a region inside the hadron, where confinement arises. From this theory, confinement 
may be considered as an intrinsic propriety of the space-time, that when interactions with large coupling constants 
are considered, generates trapping surfaces. On the other hand, no collapse for r = is expected, the uncertainty 
principle forbids it. £ as defined in (|39|l is a function of a/m, so, in Nature, a and m are such that the confinement 
conditions are filled and in a region with the observed size (some examples are shown in Table II) . 

Thinking matter as composed of small black holes may seem a weird idea, but no one has actually seen a quark, or 
an element of strongly interacting matter, and in this sense, a black hole is quite reasonable. At astronomical level 
also, it is possible to imagine systems that cannot be seen, due to the curvature caused by electromagnetic or strong 
forces, and maybe giving an important contribution to the mass of the universe. 

The most important feature of the theory, is the fact that the insertion of general relativistic aspects in the quantum 
theory generates results that are in accord with the phenomenology of the considered systems. 
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